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Abstract. We calculate the neutrino energy emission 
rate due to singlet-state pairing of protons in the neutron 
star cores taking into account the relativistic correction to 
the non-relativistic rate. The non-relativistic rate is nu- 
merically small, and the relativistic correction appears to 
be about 10-50 times larger. It plays thus the leading 
role, reducing great difference between the neutrino emis- 
sions due to pairing of protons and neutrons. The results 
are important for simulations of neutron star cooling. 

Various mechanisms of neutrino emission in the cores 
of neutron stars (NSs) arc important for NS cooling 
(e.g., Pethick 1992). In this Letter, we analyze a spe- 
cific neutrino reaction associated with superfluidity of 
protons. The reaction consists in neutrino-pair emission 
p — ► p+is+9 by protons whose dispersion relation contains 
superfluid energy gap. In theoretical studies, this process 
is treated as annihilation of quasiprotons p into neutrino 
pairs, p+p — > v+D, and is referred to as neutrino emission 
due to Cooper pair formation. Similar but more powerful 
process takes place due to pairing of neutrons. 



It is widely accepted (e.g., Takatsuka & Tamagaki 1993 



and references therein), that nucleon superfluidity in the 
NS cores occurs mainly due to triplet-state pairing of neu- 
trons and singlet-state pairing of protons. The neutrino 
emission due to singlet-state pairing of neutrons (as the 
s imple st example) has been calculated b y Flo wers et al. 
( 1976 ) and Voskresensky & Senatorov ( 1986 ) while the 
cases of triplet-state pairing of neutrons and singlet-state 
pairing of protons have been studied by Yakovlev et al. 
( |1999| , hereafter Paper I). The importance of neutrino 
emission produced by nucleon pairing has been demon- 
strated by simulations of NS cooling (e.g., Schaab et al. 



1997, Page 1998, Paper I) 



So far the neutrino energy emission rate (emissivity) 
due to pairing of nucleons has been calculated in the 
approximation of non-relativistic nucleons. We will show 
that the emission due to singlet-state pairing of protons is 
actually determined by the relativistic correction. 
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The process in question goes through neutral elec- 
troweak currents and is accompanied by emission of neu- 
trinos of all flavors. As shown in Paper I, the non- 
relativistic term in the neutrino emissivity is produced by 
vector electroweak proton current and can be written as 
(h = c = fes = 1) 



Qo 



4G F m*p F 



15tt 5 



T 7 Af„a F(y), 



(1) 



where Gf is the Fermi weak interaction constant, T is tem- 
perature, to* is an effective proton mass in dense matter 
(determined by the proton density of states near the Fermi 
level), pf is the proton Fermi momentum, and M v = 3 is 
the number of neutrino flavors. The function 



F(y) = v 2 



(e* + l) 2 ' 



(2) 



depends on y = 
gap; x = rj/T, z 



A/T, where A is the proton superfluid 

= E/T 



V^ 2 



y 2 , E is the quasiproton 



energy with respect to the Fermi level, n = vf(p — Pf), 
vp = p-p/rrip is the proton Fermi velocity. Finally, ao = 
(1-4 sin 2 9 W ) 2 w 0.0064, and W is the Weinberg angle 
(sin 2 9w = 0.23). Numerical smallness of ao for protons 
comes from their quark structure. Accordingly, the neu- 
trino emission due to proton pairing is greatly reduced 
as compared with the emission due to neutron pairing; 
for instance, one has ao = 1, for singlet-state pairing of 
neutrons. 

Let us calculate the relativistic correction Q\ to Qo 
determined by the axial-vector electroweak proton cur- 
rents. The corresponding interaction Hamiltonian is Ha = 
— Gfca J a l a /(2V2), where J a — ijj p ~f a ^ 5 ip p is the axial- 
vector proton 4-current (a —0, ...3), l a — (V^7 Q (1 + 
7 5 )Vv) is the neutrino 4-current; ip p is the bispinor wave- 
function of quasiprotons (normalized with respect to one 
particle in unit volume), ip u and %j)' v are the standard 
bispinor amplitudes of emitted neutrinos; upper bar de- 
notes a Dirac conjugate, 7° and 7 s are Dirac matrices, and 
ca is the axial-vector constant (for protons, ca — ~9A = 
— 1.26). The bispinor ip p can be presented as combination 
of an upper spinor (j> p and lower spinor Xp- Writing the 
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4-vector J a as (J°, J), we have J° = — <t> p Xp~ X p 4>p> an d 
J = —<ftp(T(f> p — Xp (Tx p , where a is the Pauli matrix. 

Let i& be the non-relativistic second-quantized Bogoli- 
ubov spinor wave function of quasiprotons in superfluid 
matter (Lifshitz & Pitaevskii 1980). For our purpose, it is 
sufficient to set cf) p = 9, Xp = — wV$/(2m p ), where m p 
is the bare proton mass. Since the lower spinor Xp is small 
as compared to <p p we can neglect the term Xp^Xp m J- 
Further considerations are similar to those in Paper I, and 
we omit the details. 

The squared matrix elements of the interaction Hamil- 
tonian summed over quasiproton spin states contain the 
tensor components I af3 = £s P i„s ( B \ ja \ A ) ( B \ J P \ A )\ 
where \A) and \B) denote, respectively, initial and final 
states of the quasiproton system (e.g., Paper I). Direct 
calculation shows that the tensor I a P is diagonal, with 
two nontrivial matrix elements: 



I t ee I m = {u'v + v'uf\p-p'\ 2 /(2ml), 

I s EE I 11 =I 22 =I 33 = 2(U'V-V / U) 2 , 



(3) 



where u = ^(E + rj)/ (2E) and v = y/(E - rj)/{2E) are 
the coefficients of the Bogoliubov transformation for an 
annihilating quasiproton with momentum p and energy 
E = y/ rf + A 2 , while vl and v' are the same coefficients 
for a second annihilating quasiproton with momentum p' 
and energy E' . The general expression for the neutrino 
emissivity Q\ due to axial- vector proton interaction Ha 
is readily given by Eq. (13) or (16) of Paper I, where cy 
has to be formally replaced by ca- Taking into account 
diagonality of I a P we obtain 
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dp dp' f(E)f(E')u> 



[k 2 I t + (3oj 2 - 2k 2 ) I s ] . 



(4) 



In this case, lu and k are, respectively, energy and mo- 
mentum of a neutrino pair (u> = E + E',k = p + p'), 
f(E) = (e B / T + 1) _1 is the Fermi-Dirac distribution; in- 
tegration over p and p' is done over the kinematically 
allowed domain lj 2 > k 2 . 

Since the proton Fermi liquid in NS cores is strongly 
degenerate, we set p = p' = in all smooth functions un- 
der the integral. The presence of energy gaps opens the re- 
action kinematically in a small region of momentum space 
where p is almost antiparallel to p'. This enables us to set 
p' = — p in all smooth functions in the integrand. Then 
Eq. (|) gives I t = 2(p F /m p ) 2 (A/E) 2 while I s -> for 
p — > p'. Expanding / s in the vicinity of p — p' we obtain 
I s = VpA 2 (p — p') 2 /(2E 4 ). After some transformations 
described in Paper I we come to Eq. (20) of Paper I in 
which again cy should be replaced by ca- Subsequent in- 
tegration over p and p' is quite analogous, and finally we 
come to the expression for Q\ which can be obtained from 
(0) by replacing ao — > Oi, where 



«i 



c 2 A v 2 [(m* p /m p ) 2 + 11/42] 



As seen from the derivation, both non-trivial tensor com- 
ponents, It and I s , contribute to Qi, producing the first 
and second terms in square brackets in Eq. (^), respec- 
tively. The relativistic correction to the neutrino emissiv- 
ity due to singlet-state pairing of neutrons has been cal- 
culated by Flowers et al. (1976). It should be the same 



for neutrons and protons since the squared axial-vector 
constant c 2 A — g 2 A is the same. However, the expression 



obtained by Flowers et al. ( 1976J) contains only the sec- 
ond term in square brackets in (|5|) . Clearly the first term 
was overlooked. 

Introducing Q = Qo+Qi an d a = ao+ii we obtain the 
total neutrino emissivity (in the standard physical units) 



Q 



4Glm*p F 
157r 5 ^ 10 c 6 

* \ 2 

mz 



{k B T) 7 Af v aF(y) 
— TlN v aF(y) 



- 1.170 x 10^ 

erg 



0.0064+ 1.60 (v F /c) 2 (ml/mp) +11/42 



cm° s 

2 



(6) 
(7) 



Let us remind a useful relationship v-p/c — 0.353 
(m p /m*) (rip/no) 1 ^ 3 , where n p is the proton number den- 
sity and no = 0.16 fm -3 is saturation nuclear matter den- 
sity. 

One can see that the relativistic correction greatly ex- 
ceeds the non-relativistic term under typical conditions 
prevailing in the NS cores. Notice also that Eqs. (0) and 
(0) neglect the relativistic corrections to Q due to vector 
proton currents. [] This is justified since the main term, 
in our case, is numerically small and unimportant. In the 
cases of any neutron pairing or triplet-state proton pairing 
the non-relativistic term is not numerically small and all 
relativistic corrections can be neglected, in practice. 

The fit expression of F(y) is (Paper I) 



F(y) = (0.602 y 2 + 0.5942 y 4 + 0.288 y 6 



x 



0.5547 + x /(0.4453) 2 + 0.01130 y- 



1/2 



x exp (-v^y 2 + (2.245) 2 + 2.245J ; (8) 

y = Vl" 7 7 (+456-0.157 t~ 1/2 + 1.764t~ iX 

where r = T/T cp < 1, and T cp is the superfluid critical 
temperature of protons. Equations (||) - (||) enable one to 
calculate the neutrino emissivity Q. 

Let us discuss the importance of neutrino emission due 
to Cooper pairing of protons in superfluid NS cores. For 
illustration, we will adopt a moderately stiff equation of 



state proposed by Prakash et al. (1988) (the same version 



(5) 



as used in Paper I) . We assume that dense matter consists 
of neutrons, protons and electrons (no muons and hyper- 
ons) . The equation of state allows the direct Urea process 

We have shown that inclusion of these corrections leads 
formally to Eq. (1), where F is given by Eq. (2) but with addi- 
tional factor C under integral: (7=1 + VpD/z 2 , D=—5zy 2 /21 
+ (10z/63)(y 2 - 3x 2 )/(e z + 1) + (10/21)( Z 2 + y 2 ) - {zm;/m p ) 2 . 
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Fig. 1. Temperature dependence of neutrino emissivities 
in npe matter at p = 2pQ. The reactions are: modified 
Urea (MU, sum of n and p branches), nucleon-nuclcon 
bremsstrahlung (NN, sum of nn, np and pp branches), 
Cooper pairing of neutrons (CPn) and protons (CPp). 
The emissivity from the latter reaction is shown twice, 
in the non-relativistic approximation (dots), and with the 
relativistic correction included (solid line). 



to operate at p > 4.64 p , where po — 2.8 x 10 14 g cm -3 is 
the saturation nuclear matter density. Consider two den- 
sities, p = 2 po and p = 5 po, as examples of the standard 
neutrino cooling (direct Urea is forbidden) and the cooling 
enhanced by the direct Urea process, respectively. The nu- 
clcon effective masses will be set equal to 0.7 of the masses 
of bare particles. Neutrons and protons are supposed to 
form superfluids (triplet-state neutron and singlet-state 
proton pairing) with critical temperatures T cn and T cp , 
respectively. 

We compare the neutrino emissivity due to pairing of 
protons with those resulting from the main neutrino re- 
actions in superfluid NS cores: direct and modified Urea 
reactions, nucleon-nucleon bremsstrahlung (as described 
in Yakovlev & Levenfish 1995 and Levenfish & Yakovlev 



1996), and Cooper pairing of neutrons (Paper I). Our anal- 



ysis shows that the emission due to pairing of protons is 
especially important for T T cp <C T cn . In this case, 
strong neutron superfluidity greatly suppresses almost all 
other neutrino reactions and Cooper pairing of protons 
becomes dominant. 

Figure 1 shows temperature dependence of the neu- 
trino emissivities at p = 2po, T cn = 6 x 10 9 K and 
T cp — 8 x 10 s K. For T > T cn matter is nonsuperfluid, 
and the modified Urea reaction dominates. For lower T, 
neutrino emission due to neutron pairing becomes most 
important but it dies out at T < 0.1 T cn . The emission 
due to proton pairing dominates aX T < T cp . Including 



P=5 Po / 
T„ = 6 x 10 9 K ^>/ / 
T rn =8 xl0 8 K / 
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Fig. 2. Same as in Fig. 1 but for p = 5 pQ. Additionally, 
direct Urea process (DU) is allowed. 

the relativistic correction increases its efficiency by a fac- 
tor of about 10. 

Figure 2 is the same as Fig. 1 but for p = 5 pq. Now 
the powerful direct Urea process is allowed and dominates 
at T ^ 3 x 10 8 K. However, in spite of its high efficiency 
in nonsuperfluid matter, it is strongly suppressed at lower 
T, so that the neutrino emission due to proton pairing 
becomes dominant. The relativistic correction enhances 
this process stronger than at p — 2p n , by a factor ~ 30. 

Concluding, the relativistic correction greatly en- 
hances the neutrino emission due to proton pairing and 
brings it closer to the neutrino production due to neutron 
pairing. This is important for NS cooling. 
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